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Abstract 

We investigate black hole formation by a spherically collapsing thin shell of 
matter in AdS space. This process has been suggested to have a holographic in- 
terpretation as thermalization of the CFT on the boundary of the AdS space. The 
AdS/CFT duality relates the shell in the bulk to an off-equilibrium state of the 
boundary theory which evolves towards a thermal equilibrium when the shell col- 
lapses to a black hole. We use 2-point functions to obtain information about the 
spectrum of excitations in the off-equilibrium state, and discuss how it character- 
izes the approach towards thermal equilibrium. The full holographic interpretation 
of the gravitational collapse would require a kinetic theory of the CFT at strong 
coupling. We speculate that the kinetic equations should be interpreted as a holo- 
graphic dual of the equation of motion of the collapsing shell. 
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1 Introduction 



Following the discovery of the AdS/CFT correspondence |J a large amount of work has 
been devoted to clarifying and extending its implications. This work has been recently 
reviewed in [[J where an extensive list of references can be found. In this context, we have 
investigated f|, Q the relation between physics in the bulk of AdS and its holographic 
dual. We have mainly focused on the formulation in Minkowski signature || [7] since 
this allows the study of dynamical processes. In particular we have been interested in the 
holographic interpretation of black hole formation in the bulk since the full understanding 
of this process is an important step towards the ultimate goal of finding a holographic 
interpretation of the black hole information problem |§ , and (hopefully) its solution. We 
also believe that by investigating the formation of black holes one will gain a deeper 
understanding of general features of holography. The process that we are focusing on is 
the spherical collapse of a thin shell of unspecified matter in the bulk. A special case, the 
spherical collapse of a dust shell in AdS3, was analyzed by Peleg and Steif in ||[]. In the 
holographic dual, a flat space field theory with no gravity, the bulk process is expected to 
correspond to the time evolution and thermalization of an initial out of equilibrium state 
]7|, [□], [L2|, [13| . The collapse of the shell due to gravity describes the approach towards 
equilibrium, while the final formation of the black hole corresponds to thermalization. It 
is intriguing how a gravitational process in the bulk is encoded in a kinetic process in 
the holographic dual. This, we believe, suggests a deep connection through holography 
between gravity and kinetic theory of the boundary theory. 

In [|J we developed Green function techniques that allowed us to obtain the holographic 
"images" of various objects in the bulk. In |5j we used these techniques to begin our study 
of black hole formation. As a first step we studied the case of a very large shell of matter 
which is just beginning to collapse. In this case we could neglect the motion of the shell 
and use a quasistatic approximation. The motivation for the study was to find out how 
the scale-radius duality f7|, [12|, [U| encodes the radial size of a spherically symmetric 
shell, and to find properties of the off-equlibrium state in the holographic dual. We found 
that the presence of the shell manifested itself in the form of a series of poles in the 
propagator. The location of these poles depends on the radius of the shell, in agreement 
with the scale-radius duality. A slight surprise was that the poles also had an imaginary 
part. Hence they can be interpreted as unstable collective excitations, "shellons" (since 
they exist because of the shell), in the off-equlibrium state. 

In this paper we will carry the analysis one step further. While we are not yet able to 
address the time-dependence of the collapse, there is another regime where a quasistatic 



4 The other extreme case, collapse of a solid disk of dust in AdS3 was investigated by Mann and Ross 



1 



approach can again be adopted: the very last stage when the shell is approaching the 
horizon radius rjj of the black hole about to form. There are two ways to justify the 
quasistatic approximation. 

First, in the final stage of the gravitational collapse the shell appears to "freeze" in 
the frame of an asymptotic observer. The approach of the radius of the shell r s to the 
radius of the horizon r# is exponentially slow, 

where the characteristic timescale th is inversely proportional to the Hawking tempera- 
ture[], 

Therefore, the quasistatic approximation is applicable for evaluating the propagator as 
long as we focus on energies which are larger than the inverse characteristic timescale, 

\u\ > — ~ T H . (3) 

At the same time, this is the region of interest for investigating the approach to ther- 
mal equilibrium. At equilibrium, the propagator becomes thermal, with a characteristic 
infinite sequence of thermal poles 

u n — i4irnTH + const. , n = ±1, ±2, . . . . (4) 

When the radius of the shell is slightly larger than the horizon radius, the boundary theory 
is slightly out of equilibrium, and one finds a correction factor to the thermal propagator 
which has additional (shellon) poles and zeroes. They characterize the deviation from 
the equilibrium; the poles and zeroes flow together and disappear during the approach to 
equilibrium. 

Second, one can avoid the restriction to large frequencies by investigating a slightly 
different way of approaching the equilibrium. In the bulk one can consider a sequence 
of static spherical shells with a constant total ADM mass, but with progressively smaller 
radii approaching the horizon radius. Instead of starting with a shell at rest and letting 
it collapse freely, this corresponds to slowing down the collapse by applying an external 
force so that the process becomes quasistatic, while increasing the amount of matter in 
the shell. The ADM energy needs to be kept fixed to ensure that the shell approaches its 
horizon radius, therefore the loss of kinetic energy needs to be compensated by adding 
matter to the shell during the process. In this case one obtains more detailed information 

5 This is true for the production of a generic black hole in a generic dimension. 
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about the propagator in the boundary theory. Presumably, the quasistatic approximation 
amounts in the boundary theory to replace the kinetic theory by a thermodynamical 
treatment of small deviations from equilibrium. In thermodynamics the evolution of a 
system is always considered to be quasistatic, and small deviations from equilibrium can 
be arranged to be such by applying external forces. The evolution is then determined by 
the gradient of entropy depending on a small parameter characterizing the deviation from 
equilibrium, which should correspond to the applied force balancing the gravitational force 
in the bulk. In contrast, recently |TJ| the equations of motion in the bulk were related to 
renormalization group equations in the boundary. This is a good interpretation for the 
static background geometry but general relativity also determines the time evolution of an 
object in the bulk - this is a dynamical process related to the dynamics of the boundary 
as discussed above. 

The paper is organized as follows. In section 2 we review our previous work and 



compare our results with other recent work on black hole formation [17], |18, 19, ^C§ . The 



analysis of the last stages of black hole formation in the case of AdS3 is given in section 
3. Section 4 contains a discussion of results and some ideas for the future. 



2 Derivation of the correlators 



We begin by briefly reviewing (following [21, 22] and as described in ||) the essential 
steps of the computation of a two point function of a boundary operator coupling to a 
bulk scalar field in an asymptotically AdS spacetime. A scalar field in the bulk satisfies 
the classical equation of motion 



'-9 



<f>{t,x,r) 







(5) 



where g^ v is the bulk metric and the mass term m 2 contains the mass of the field, the 
coupling to the curvature scalar^], and contributions from a Kaluza-Klein reduction on a 
compact internal space, if the original spacetime is higher dimensional. At first, we choose 
the desired coordinate system for (a part of) AdS space. We will work in Minkowski signa- 
ture. Depending on the choice of coordinates, we then impose the appropriate boundary 
condition for the bulk field in the interior. For example, in AdS with global coordinates 
one requires that vanishes at the origin. In Poincare coordinates or in the case of a 
black hole, the appropiate boundary condition at the horizon is an ingoing wave when 
the imaginary part of the energy is positive and an outgoing when it is negative. The 

6 When we consider e.g. a spacetime of a collapsing shell in AdS,j+i with d > 2, the exterior metric 
has a nonconstant scalar curvature. In such cases we assume that the scalar field is minimally coupled 
so that m 2 remains a constant. 
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propagator then has a cut on the real axis where the imaginary part changes sign. Let 
t, x denote the coordinates that parametrize the boundary, and let r denote the radial 
coordinate. We perform a Fourier transformation to momentum space in the boundary 
coordinates, so the interior solution is (f>(u>, k, r). This solution is uniquely determined (up 
to an overall normalization) by the boundary conditions discussed above. 

We then study its behavior near the boundary. In Minkowski signature, there are two 
kinds of bulk field modes, normalizable and non-normalizable[]. In momentum space, we 
denote the former by <j>^(u, k, r) and the latter by (f)(~\uj, k, r). We normalize the modes 
such that their asymptotic behavior near the boundary is 

(±) (cj,£,r) = r~ A± • 1 (r^oo) (6) 

where 

d 



A± = -(d ± Vd 2 + Am 2 ) = - ± u . (7) 

Now, in the asymptotic region, the interior solution 4>(u, k,r) becomes a linear combina- 
tion of the normalizable and non-normalizable solutions, with an asymptotic behavior 

x, r) « r~ A+ C(u, k) + r~ A -D(uj, k) . (8) 

The coefficient D(u,k) is the Fourier transform of the boundary data 4>o(t,x) which acts 



as a source term for a dimension A + operator 0(t,x) in the boundary theory EI], |22 
Next, we rewrite (H) in the form 



<P(uj, k, r) w [r- A+ G(u;, k) + r~ A -\ 4> (uj, k) , (9) 

where 

G(w ,£) = ^|. (io) 

D{u>, k) 

One can then check that the (finite part) of the two point function of the operator O at 
the boundary has the momentum space formP] 

(0(u, k)0(u', k')) = ( + ~ ~ \ 5{u + J)8{k + k')G{u, k) . (11) 

In the case that the geometry has a horizon, the propagator, as a function of the 
frequency u, has a cut on the real axis which arises from the boundary conditions as we 
already discussed. In that case one can analytically continue the propagator from the 



7 The Minkowski signature version of AdS/CFT duality was first investigated in j(| |7j. 
8 The evaluation of the overall coefficient involves a subtlety, see ^3|. In the remainder of the paper 
we will suppress overall coefficients as they do not affect the pole structure of the propagator. 
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upper half-plane into the lower half-plane. This analytic continuation can have poles on 
the lower half plane signaling the presence of unstable quasi particles. 

Next, we compute the propagator in the presence of a thin matter shell of radius r s 
in AdSd+i space. We are going to restrict attention to a quasistatic approximation: we 
assume that the shell is collapsing very slowly (as measured in the asymptotic time), 
so that it appears to be static^. The static case has an interest also of its own since it 
provides a series of states which are closer and closer to equilibrium. Understanding how 
these states differ from the black hole from the boundary point of view, gives an insight 
into the boundary processes at work. Furthermore, the usual thermodynamical approach 
to out of equilibrium states is to consider states which are at equilibrium but for a variable 
that is externally kept fixed. In this case the variable is the radius of the shell. 

In the interior of the shell, the metric is the AdS metric in global coordinates, and the 
exterior metric is the AdS black hole metric. Both metrics can be expressed in the form 



dv 2 

ds 2 = -f{r)dt 2 + — - + r 2 ^ d _! , (12) 
fi r ) 



with the function f(r): 



f( r )-i fi(r) = l+r 2 if r < r s 

J[r) \ /2 ( r ) = l--^ + r 2 if r>r s , [i6) 

where the parameter \x is related to the ADM mass of the black hole. The radius of the 
horizon of the black hole, r#, is determined by solving the equation fiiru) = 0. (We 
are using units where t, r and the angular coordinates are dimensionless and the AdS 
radius is set to R — 1 . The dimensions can be restored by replacing everywhere r —>■ r/R 
etc. With the metric (|T2"D, the scalar field equation (|5p can be reduced (after a Fourier 
transformation to momentum space in the boundary coordinates) to the radial equation 

-^d r (r^ 1 / W) <f>(u, k, r) + (j^ - ^ - m 2 ) 0(c, k, r) = , (14) 

where k 2 = 1(1 + d — 2), I — 0, 1, 2, ... is the eigenvalue of the Laplacian on the (d — 1)- 
sphere. The interior and exterior solutions <pi and <p2 of the radial equation must be 
matched across the thin shell using flux conservation. The matching conditions are 

4>l\r=r a — 4>2\r=r a (15) 
fl(r)d r 4>i\ r=rs = f2(r)d r 4> 2 \r=r a ■ (16) 



9 Of course, a complete treatment should take into account what happens when the shell is collapsing 
rapidly. This can be expected to give rise to particle production in the boundary theory even if we neglect 
the particle production in the bulk. We will leave these more demanding issues aside here. 
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To find the propagator G(u, k) in the boundary, we first expand the exterior solution in 
normalizable and non-normalizable modes, 



cj) 2 (uj, k, r) « r-*+C(uj, x) + r^-D{u, k) , (17) 

then substitute (|T7| ) in the matching conditions above, solve for the ratio C/D in terms 
of (f>i and (f) 2 , and substitute into equation ( |T0"D which gives the propagator. The result is 



Vf2d r (p2 - n fldr 



In the above, all quantities are evaluated at the radius of the shell r s . The analytic 
structure of the resulting propagator gives information about the spectrum of excitations 
in the boundary theory. In particular, we are interested in the poles of the propagator at 
zero transverse momentum (k = 0) in the complex cj-plane. Poles at positive real cu-axis 
correspond to masses of stable composite objects created and annihilated by O. However, 
we will typically find complex poles at 

w n = M n - l -T n (19) 

corresponding to unstable resonances created and annihilated by O. In a previous paper ||, 
we used the name "shellons" for the resonances. The real part M n gives the mass of the 
resonance, and the imaginary part T n is the width which is the inverse of the lifetime. 
There can also be poles in the imaginary cu-axis, for example when the black hole has 
fully formed and the boundary theory has thermalized. Then the boundary propaga- 
tor becomes a thermal propagator, and the imaginary poles are related to periodicity in 
imaginary time. 

From , we see that the poles typically correspond to the zeroes of the denominator 
(unless they cancel against those of the numerator). The vanishing of the denominator 
means that the interior solution <fri matches completely with a normalizable exterior so- 
lution (f>2 ■ Thus, the total mode is normalizable - unless u is complex, in which case it 
will be quasinormalizable. For a better understanding of the situation, it is instructive to 
map the radial equation into the form of a Schrodinger equation and consider the resulting 
quantum mechanical analogue. We rescale the field (p as 

<f>(u, k, r) = A- d+1 V 2 ifj{uj, k, r) (20) 

and use a "tortoise" coordinate 

W/M' (21) 
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Note that we have to use fi in the interior and fa in the exterior, so we must add an 
integration constant to ensure that r* is continuous across the shell. The range of r* is 
finite. For example, if d = 2, we get 



2r H \r-r H J 



T > Ts (22) 
arctan(r) — r < r < r s 



where ^ 

r = arctan(r s ) + - — In ( — — — ) (23) 
2r H \r s -r H J 

and the full range of r* is the finite interval — r < r* < 0. In the following we will mostly 
be considering the case where the radius of the final black hole, r#, and therefore also the 
radius of the shell, r s , are very large compared to the AdS radius (r s > r# ^> 1 in our 
units). The radial equation in the interior will then reduce to that in a Poincare patch of 
AdS space (see f5|), and the two equations above become 



2r H \r-r H J 



where 



(24) 

L — - — fQ < r < r s 

(25) 

2r H \r s -r H J r s 
The full range of r* is now the infinite interval — oo < < 0, corresponding to the infinite 
AdS throat. 

With the above rescalings, the radial equation (13) reduces to the form of a time 
independent Schrodinger equation, 



- d 2 r V + V(n)ip = u; 2 ij (26) 



with a potential 



= (d-l)(d-3)p + (d-l)d r J + (fc 2 + m 2 r)/ ^ 
4r 2 2r r 



where r = r(r*) from inverting (^T ) and / = f(r(r*)). Thus the problem of solving the 
radial equation becomes analogous to a one dimensional quantum mechanical scattering 
problem. In Figure 1 a typical profile of the potential is depicted. It blows up at the 
boundary and has a small barrier at the location of the shelQ There is also an infinite 
well at = — r , but when we consider a large black hole the well is pushed down the 

10 If we were relax the quasistatic approximation and allow for the time dependence of the collapsing 
shell, we would need to consider quantum mechanical scattering with a time dependent potential: the 
barrier would be moving to the left. 
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Figure 1: Potential in the presence of a shell. 



AdS throat towards r* = — oo and plays no role in the calculation. Now let us return to 
the poles of (0). As we discussed above, the appropriate interior boundary condition for 
0i is that it is a purely ingoing wave. On the other hand, the normalizable mode in the 
exterior is a superposition of an outgoing wave and an ingoing wave, due to the reflection 
from the boundary of AdS. For example, it is simple to check this by looking at the high 
frequency limit of the normalizable mode (p^ 0, 0, which looks like 



,M+D/2 COS (^ + O ) 



(2c 



where 9q is a 
superposition 
of the barrier 
is a purely inj 
impossible to 
complex energ 
frequencies u>. 
( HH) takes an 



constant phase shift. For the "wavefunction" ip this means that it is a 
of ingoing and outgoing waves with equal amplitude on the right hand side 
corresponding to the shell (see Figure 1) and on the left hand side there 
;oing wave. Thus there is total reflection together with transmission - it is 
satisfy both conditions with real energy levels. However, it is possible with 
;y levels u 2 . This is another way to see why the poles appear at complex 
Incidentally, note that using the tortoise coordinate fl2"ID , the propagator 
even simpler form - it is just a ratio of two Wronskians 



G{u,k) 



(29) 



evaluated at r* = r*(r s ). 
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Let us compare our approach with other recent work. Horowitz and Hubeny fl9fl 
investigated the quasinormal modes in the background of a black hole0. They study the 
equation (JT4]) with the black hole metric (/(r) = /2(?"))- The boundary conditions in 
|19| for a scalar field are very similar to ours: one takes the field to be a superposition 



of outgoing and ingoing modes near the boundary (normalizable modes), whereas the 
interior boundary condition in is that there should be only ingoing modes near the 
horizon of the black hole. In our case with a shell, the interior metric is that of a pure 
AdS space. However, if the radius of the shell is larger than the AdS space (r s 3> 1 in the 
present units), a Poincare horizon appears at r — > (r* = — oo) due to the approximation, 
and there is a natural choice between ingoing or outgoing modes, which is reflected in a 
cut appearing in the propagator^ of the boundary operator O. Thus, on a formal level 



the modes we discussed in || are analogous to the quasinormal modes of black holes |19 



The complex spectrum arises because of the choice of boundary conditions, and the main 
difference between [19] and || is that the horizon of a black hole has been replaced by a 
Poincare horizon. However, the physical interpretation appears to be somewhat different. 
The horizon of a black hole is a real boundary through which matter can fall in and be 
lost, while the appearance of a Poincare horizon in the above is merely an artifact of the 
approximation since we started from global AdS space. In reality we would expect an 
ingoing wave to be reflected back from the origin of the AdS space and not to be lost. 
What happens is that if the radius of the shell is much larger than the AdS radius, r s ^> 1, 
and the frequency is complex (with ^slj < 0) the reflection is exponentially attenuated, 
so that it can be suppressed. There are two other ways to arrange for the ingoing mode 
from the shell not to come back at all. One way is simply to consider a shell falling into 
an existing (small) black hole instead of creating the black hole in the collapse. Another 
way is to consider late times in the collapse when the shell is approaching the radius of 
the horizon of the black hole about to be created. From the space-time diagram of black 
hole formation in a spherical collapse, one can see that the horizon forms even before the 
black hole has formed. Therefore, the ingoing modes from the shell will be lost into the 
region of trapped surfaces and never come back. 

Another related work is by Giddings and Ross |1| who investigated a collapsing shell 
of D-branes. In their case, contrary to ours with an unspecified matter shell in AdS space, 
the geometry in the interior of the D-brane shell is flat. Hence they do not have an infinite 
throat, where objects falling in could get lost. Instead, their space is cut off at a finite r*. 
As a consequence, when they investigate the poles of the propagator, they initially obtain 
a discrete spectrum of stable excitations with real energies, much like the global modes in 

11 Related work can be found in p4| . 

12 The choice between the two sides of the cut corresponds to a choice between a retarded or advanced 
propagator. 
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our case. However, the shell itself can absorb incoming waves, and after considering this 
effect they also obtain imaginary contributions to the energy spectrum. So the absorption 
by the shell plays a similar role as the absorption by the infinite throat in AdS space in 
our case. 



Furthermore, Balasubramanian and Ross JT7j have investigated black hole formation 
in AdSs by point particles. The bulk solution was investigated by Matschull ||25|| , and 
Balasubramanian and Ross investigate the use of 2-point functions in the boundary theory 
to keep track of the positions of the point particles as they approach each other. In this 
case one does not have the complication of spherical symmetry, the point particle positions 
can be read off from kinks in the 2-point functions in coordinate space. 

Finally, we would like to mention a slightly different line of work on thermalization in 
the context of black holes and AdS/CFT correspondence. Kiritsis and Taylor [26] investi- 
gated D-brane probes falling into a black hole. In the process, the potential energy of the 
probe is converted into heat, which then is absorbed by the black hole. By investigating 
the bulk action of the probe, one can e.g. derive the equation of state for the black hole. 
But, on the other hand, one can study the dual gauge theory effective action of the probe 
at finite temperature and try to reproduce predictions from bulk calculations; this was 



investigated in 27 



3 Black hole formation in AdS3 

In this section we examine the boundary propagator in the case that a shell is quasistat- 
ically forming a black hole in AdS3 space. As in ||, we assume that the resulting black 
hole has a much larger radius than that of the AdS space: r# ^> 1 in the present units. 
The problem is to find the modes 0i,02 which are needed in the propagator formula 
fll8|) . To find the interior mode, note that in the region r > r s > 1 we can approximate 
/i(r) = 1 +r 2 « r 2 . Then the radial equation becomes the mode equation in Poincare co- 
ordinates, and the mode solutions will be Bessel functions <fi ~ (l/r)J± u (yuP — k 2 /r). As 
discussed in section 2, in Poincare coordinates the natural choice for the interior boundary 
condition is between ingoing and outgoing modes, and we choose the former, as we did 
in pi. However, in H we made an additional approximation and focused on the large 
frequency limit u/r ^> 1. Here we will not make that approximation, so the infalling 
mode is the following linear combination 

tl = Bl l l jJ^EE).^j(^EE) ) , ,30) 

r \ r J \ r J 

where v = a/1 + m 2 and c% is a normalization constant (which will drop out in the answer 
for the propagator). 
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f ] = (u - lTu-^-«F{a + ^± a + ^; 1 + u; u' 1 ) , (31) 



Then, we need the exterior modes 02 in BTZ coordinates. In our previous paper ||, 
we used a WKB approximation in the large frequency limit to obtain the modes. Here, 
we will instead use the exact solutions of the wave equation p8| , the normalizable and 
nonnormalizable modes can be found in ||. If the shell is large (r s ^> r ), we use 

A± A 

where u = r 2 /r 2 Hl a = iu;/(2r^) and A± were defined in Substituting (|30|) and (0) 
into ([18|) , we obtain the result for the propagator of the operator O on the boundary in 
the presence of a large matter shell in the bulk. Using the latter result we checked that 
we reobtain our previous approximate result for the poles in the r s 3> ru, u ^> r s limit: 

3 v , ir s , / Afmr 2 \ , , 

«» = Tr.(n+- + -)--^ln^— . (32) 

This is sufficient to illustrate the main point: the presence of the spherically symmetric 
shell in the bulk is encoded in a tower of resonances in the boundary theory, with the 
above spectrum from which one can read their masses and decay widths. An intriguing 
thing about the shellon poles becomes apparent if we focus on the UV limit |u;| ^> 1. 
By the scale-radius duality one might have thought that we in the bulk are probing only 
the geometry very near the boundary. But that is a region outside of the shell, so the 
geometry is equal to the black hole geometry. On the other hand, near the boundary the 
black hole metric reduces to the same form as the Poincare metric. From these arguments 
we might have expected that the analytic structure of the propagator in the UV limit is 
either similar to that of the thermal propagator, or the Poincare propagator 13 . On the 
other hand, from (^) we know that there are poles for arbitrary large uj arising from 
the shell. A possible interpretation is that the poles of the propagator do not only carry 
local information about the geometry of spacetime near the boundary - they also carry 
information about an event that happens deep in the bulk inside the shell: the absorption 
of the ingoing wave into the (fictitious) horizon, which was used as a boundary condition 
in the evaluation of the propagator. 

Rather than present more details and more accurate results for the poles in the large 
shell case, we turn our attention to the other extreme of the collapse of the shell: the 
approach to the horizon radius ((r — r H )/r H ^ 1) where we expect to see thermalization 
on the boundary. 

As the shell is approaching the horizon, we choose the following expressions for the 
exterior normalizable and nonnormalizable modes |§: 

# } = A±(u-l) a F(a + ^,a + ^--2a + l-l-u) 



13 



We would like to thank Per Kraus and Sandip Trivedi for pointing out this issue to us. 
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+ B ± (u-l)- a F(-a + ^ J -a + ^;-2a + l;l-u) , (33) 

where 

r(i±i/)r(-2a) 

A± " + ^) (34) 

fl, = £(1±™ . (35) 

and we used a as defined in ([II]). We then plug the interior mode fl3"0| ) and the exterior 
modes (|33| ) into the main propagator formula (0), and study its analytic structure. 
To simplify matters, we again focus on the propagator at zero transverse momentum, 
G(u,k = 0). 

Let us begin with some general observations. As the shell approaches the horizon, 
the boundary theory gets closer and closer to thermal equilibrium. Therefore, in the 
limit r s — > r#, G(u, k) should become the thermal propagator which we evaluated in our 
previous work ||. The thermal propagator (which we denote by Gt{uj, k) from now on) 
has a characteristic infinite sequence of zeroes and poles, at zero momentum they are 
located on the imaginary axis on the complex tu-plane. Before the limit r s — > r#, when 
the shell is at a small distance from the horizon, the boundary theory is slighty off thermal 
equilibrium. In that case, we find it useful to write the propagator in a f actor ized form 

G(u, k = 0) = G t (uj, k = 0) x H(cu, k = 0) , (36) 

This is useful for two reasons. First, if r s — > r# and Imu > 0, then H(u, k = 0) — > 1 
and H encodes small deviations from the thermal equilibrium when r s ^ Th- Second, 
when Imu; < 0, H(u, k = 0) has poles which are then poles of the propagator. The 
expression can be misleading if used on, for example, the imaginary axis, since Gt has 
poles at u — —HtttiTh (Th is the temperature of the black hole) whereas G does not 
since they are canceled by H. These thermal poles appear in G only after taking the limit 
r s ~~ * r H- Formally one has to consider Imu > 0, take the limit r s —>■ in which case 
H(u,k = 0) — > 1, and then analytically continue below the real axis. Since then G = Gt 
they also have the same poles. The analytic continuation is necessary since when r s — > rjj 
a new cut arises as an accumulation of poles and zeros. In fact it can be said that the cut 
due to the boundary conditions at the Poincare horizon is replaced by another one due 
to the new boundary conditions at the horizon of the black hole which forms. 

After discussing this point, we try to analyze the propagator without any approxima- 
tions for the interior and exterior modes (|30| ) and (|33|). We assume that we are initially 
on the upper half cu-plane. Then, one can see that the leading terms of the exterior modes 
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2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 
Re 

Figure 2: Flow of the first three poles of H in the complex uo plane as the shell contracts 
towards the horizon. In the numerical calculation, we used the parameter values v = 1.5, 
Th = 10 (horizon radius), and the initial radius of the shell is r s = 10.25. 



have a coefficient A±. We can pull out these coefficients in the numerator and de- 
nominator of G, whereafter they will give us the "thermal" factor Gt- In other words, 
now 

A + Vr(-^ + % 



gvum)) - - -( ;,; A ; ■} (3D 

and 




H[U,0)= 7-T 7-T . (3£ 

kf2d r (^/A + ) - {^/A^fA ' 

The "correction factor" H{u,0) is too complicated to allow us to find an analytic expres- 
sion for the poles. Therefore we have evaluated them numerically, using Maple V version 
5. Figure 2 below depicts the flow of the first couple of poles of H as the shell radius 
r s approaches the horizon radius r#. Numerical investigations tell us that the poles of 
H start out in the lower half complex cu-plane (where the poles were found as the shell 
was still large), and then move towards the real axis and the origin. Specifically, as r s is 
sufficiently close to ro, the real and imaginary parts of the poles u become much smaller 
than rx,r s . Further, the imaginary part decreases more rapidly than the real part, and 
quickly becomes many orders of magnitude smaller. This is illustrated in Table 1 below 
which lists a few characteristic values for a pole as r s approaches th- 
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1U.U4-Z.4( l 


lU.UlO 


0. ( O-U. ( Ol 


1 n nnon 

1U.UUZU 


o.4y-u.ooi 


10.00024 


4.65-0. 18i 


10.000031 


4.03-0. lOi 


10.0000020 


3.43-0.049i 


10.00000024 


3.08-0.030i 


10.000000030 


2.80-0.019i 



Table 1. Pole values uj as r s approaches th = 10. (Using z/ = 1.5.) 

We now turn to search for a simplified expression for if, and begin by making some 
approximations for the interior and exterior modes. Let's first consider the exterior modes 
4>2 • As (r s — rn)jrn becomes small, u = r 2 s /r 2 H approaches 1 and we can do a series 
expansion in (u — 1) obtaining 

0? } « {a± [(« - ir + - ir +1 ; 



+ B± [(« - iy a + - 1)- Q+1 )] } . (39) 



Similarly, we obtain a series expansion for the combination /^r^^ which also appears 
in the equation (|I8"D for the propagator, 



/ 2 9 r 0? } « 2r {A ± [(u - l) a + - l) a+1 ) 



fin 



l)- a + G(( M -l)- Q+1 )]} . 

(40) 

Recalling that a = uo/{2th) h immediately follows that (u — l) a is the leading term if 
Imtu > whereas (u— l)~ a dominates for Imu < 0. Replacing in equation ( |18| ) results in 



lim r 



G(u,r s 



t 



if Imu; > 
if Imu; < 



(41) 



Note that independently of the function 0i, i.e. independently of the boundary conditions 
at the Poincare throat, the propagator develops a cut with the right boundary conditions 
at the horizon of the black hole. At the cut the imaginary part of the propagator changes 
sign since A± = B±. 

After checking that in the limit r s — > r# the propagator reduces to the black hole 
propagator, we proceed to compute the corrections when r s is close but not equal to Th- 
It is easy to see that if the dominant term is for example (u — l) a then no corrections are 
obtained from keeping terms of order (u — l) a+n unless the term (u — l)~ a is included. 
The same occurs when Imu < and (u — l)~ a dominates. Therefore, to compute the 
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corrections we use the approximations 

# } « {A ± (u-ir + B ± (u-l)- a } (42) 
« 2ra{A ± ( M -ir- J B ± ( M -l)- Q } (43) 

in the equation for the propagator. 

With the above approximations, the expression for the correction factor H becomes 

- a](u - l) a -[■**■ + a]^(u - l)- a 

H(u, 0) « ^ - - tr» L_ (44) 

[£-a](u-l)«-[£ + a]£(ii-l)-« 

where a = 9 r ln0i. With this expression it is very easy to see that there are no poles in 
G. In the upper half plane, where the imaginary part of to is positive, a is negative and 
G — > — as r s — > r#, (w — > 1). This expression has no poles where it is valid, i.e. in the 
upper half plane. Similarly, in the lower half plane we have G —>■ — jf- which again has 
no poles where it is valid. Let us now return to H . To find the zeroes and poles of H , 
we have to solve for the values of u where the numerator or the denominator vanishes. 

We aim here to confirm what we saw numerically before, namely that the poles and 
zeros accumulate on the real axis giving rise to a cut. In that case what has to happen is 
that the poles go to zero as r s — > r# (see appendix for an example). To check this fact we 
only need to study the propagator for frequencies such that uo/r s « 1 which is fortunate 
since in that regime the equations simplify^ in two ways. First, the mode in the interior 
of the shell, 4>\ can be expanded as: 



1 f (uj/2r) u (cu/2r) v+2 

7^ jr(i/ + i) ~ i> + 2) 



_ e VKU 



[u/2r)- v (u/2r)- u+2 ^ 



r(-z/ + l) r(-i/ + 2) 



(45) 



Keeping only the leading term when u>/r s <C 1 it follows that 



" C ' 1] r~ v -^ => a = d r hi(f> 1 = i/--. (46) 



r(-z/+ 1) ' 2 

Furthermore, we are able to do one more approximation: 

B± r(2a)r 2 (-« + A±/2) 
~ r(-2a)r 2 (a + A±/2) 

« -1 + 2#(%) + l\— + ■ ■ ■ « - e - 2i[ ^ ( ^ )+7] ^ . (47) 
2 r# 



14 If w is small we are outside the validity of the quasistatic approximation for the collapsing shell case, 
but here we want to check that all poles converge towards zero. Moreover, the approximation is still valid 
if we consider the other process which was discussed in the introduction, the "slowed down" collapse at 
fixed ADM energy. 



15 



where the last exponential is introduced for convenience. Then, the two equations (for 
poles and zeros) can be written in the form 



or equivalently 



with 




(48) 



(49) 



x = u/r H , (50) 
6± = 1b(u-1) + 2[V>(^)+7] • (51) 

As can be easily seen by plotting the functions involved, the equation (|49| ) has an infinite 
number of real solutions which can be labeled by an integer n. For 1 < n <C hi(( r s — 
r 2 H )/r 2 s ) at leading order the solutions are 

(52) 



ln[^]+2[^)+ 7 ] ' 

with n an integer, corresponding to values on the real u axis. This is our result for the 
poles and zeroes of H in the limit r s — > r#. If we choose A + in (|52"D , we obtain the 
poles, if A_, we obtain the zeroes. These two do not agree as long as A + ^ A_, which 
is true for v = a/1 + m? > 0. Note that the solutions are purely real. We commented 
earlier that as r s — > r#, the poles are complex, but flow towards the real axis and the 
imaginary part quickly becomes many orders of magnitude smaller than the real part. 
What happened above is that the imaginary part was too small to show up without 
keeping more sub leading terms. In fact keeping the term of order (u/r) u in ([15|) gives 
rise to a tiny imaginary part, with a sign depending on the boundary conditions used. 
In the extreme limit, the logarithm in the denominator of (|52| ) begins to dominate more 
and more over the A-dependent term, so the poles and zeroes approach each other, finally 
cancelling out and converting into a cut on the real axis. In the appendix we give a simple 
example of how zeroes and poles can form a cut. 



4 Discussion 



In the previous sections we have shown how, as the shell moves closer to its own horizon, 
the poles and zeros of the boundary propagator in the complex u plane accumulate giving 
rise to a cut. The appearence of the cut is the signal in the boundary theory of the 
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formation of the black hole horizon in the bulk. It is related to the fact that at the horizon 
different boundary conditions (ingoing or outgoing waves) are appropriate on different 
sides of the real cu-axis. In our case these boundary conditions appear automatically 
when the shell collapses. Once the cut has formed it is meaningful to perform an analytic 
continuation of the propagator below the real axis and then new poles can appear. In 
AdSs only the thermal poles were found, but if we had performed the calculation in AdSs 
then the poles computed in [[19] would have appeared after the collapse of the shell. Let 
us keep in mind that when we talk about collapse and motion of the shell, it is only 
within the quasistatic approximation which from the boundary perspective is more akin 
to a thermodynamic treatment than a true dynamical description. 

An open question is that there appears to be two relevant timescales as the shell is 
approaching the horizon. On one hand, r, the inverse of the imaginary part of a pole of 
H, gives a timescale corresponding to the lifetime of the corresponding resonance. On the 
other hand, the equation of motion for the shell (in this region) was approximately given 
by the equation (H) with a characteristic timescale 

r H = A-kT h ~ r H . (53) 

It remains to be understood properly which one is the relevant timescale for the thermal- 
ization in this case: the lifetime of the resonances, or the time for the collapse, and how 
the latter could be found in the boundary theory. 

Our present understanding is then the following. In the beginning, as the shell is 
very large, and just begins to collapse, the boundary theory is very far from the thermal 
equilibrium. The propagator reveals this by having a pole structure which is very different 
from a thermal propagator. The poles correspond to an infinite tower of resonances in 
the boundary theory with masses proportional to the shell radius, and lifetimes inversely 
proportional to it. As the collapse begins to speed up we cannot apply the quasistatic 
analysis - there will be dynamical processes in the boundary theory which are beyond our 
control at this stage. However, as a first crude step, we can jump to the very end of the 
collapse when a far-away observer at a fixed radial distance would see the collapse to slow 
down and "freeze" as the shell is approaching the horizon radius. At that stage, as a rough 
approximation, we again apply the quasistatic analysis, in the large frequency domain. 
Now the boundary theory is approaching thermal equilibrium with deviations encoded 
in the propagator. The poles and zeroes of the propagator flow towards the real axis 
of the complex u;-plane, and finally degenerate into a cut with the propagator being the 
thermal one. The full propagator is then defined either on the upper or lower half plane, 
corresponding to a choice between a retarded or advanced propagator. A continuation 
across the cut then reveals the thermal poles on the other half-plane. I.e. G = — 4=- is 
analytically continued down through the cut into the lower half plane where one finds 
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the thermal poles. So the propagator has, as promised, become a thermal (retarded or 
advanced) propagator when the shell is infinitesimally close to the horizon and Hawking 
radiation begins to leak out and fill the anti-de Sitter space so that a thermal equilibrium 
is established. 

We emphasize that the propagator is just a measure of the properties of the state of 
the boundary theory. The analysis does not tell us anything about how the thermalization 
happens - there is no dynamical information about the boundary theory. This question is 
tied up with going beyond the quasistatic approximation, where dynamical issues become 
relevant. What one would need at that stage, is a kinetic theory for the thermalization 
of the strongly coupled boundary theory. This seems to be an untractable problem at 
the moment. One can speculate that the equations governing the thermalization of the 
boundary theory would be a kinetic theory equivalent of the equations of motion of the 
collapsing shell in the bulk. That is, if one were to understand the kinetic theory, one 
should be able to reproduce the equation of the motion of the shell. Further, it has 



been suggested before E9, 30, 31, 32| that a holographic description of even the simplest 



dynamical processes in the bulk can be very non-local - even a slighest innocuous looking 
approximation on the boundary side might completely distort the bulk interpretation. The 
same could be true about the kinetic theory / shell collapse duality. A small approximation 
to the kinetic equations might completely mess up the bulk description of the collapse. 
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Appendix 

In section 3 we discussed how the series of poles and zeros of the propagator accumulate 
on the real axis giving rise to a cut. To clarify how this happens we give in this appendix 
a simpler example in which this occurs. Consider the function 

(54) 

where A is a positive real number ans z G C. The function G(z) has poles at z — —n/X 
and zeros at z = — (n + 1/2)/ A with n a positive integer. When A — > oo the poles and 
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zeros accumulate on the negative real axis. In that limit and for — 7r < arg(z) < 7r the 
Stirling approximation gives 



Um GW = lim = lim ^ (A^e-^ = 1 



The resulting function, 1 / yfz has a cut which, as follows from the calculation, is placed on 
the negative real axis. In fact this example arises when taking the large frequency limit 
of the AdS propagator in global coordinates. In that case one recovers the propagator in 
Poincare coordinates which has a cut on the positive real axis. 
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